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Abstract 

Global Navigation Satellite Systems (GNSS), such as the Global Positioning System (GPS), are 
among the most important sensors for movement analysis. GPS is widely used to record the tra¬ 
jectories of vehicles, animals and human beings. However, all GPS movement data are affected by 
both measurement and interpolation error. In this article we show that measurement error causes 
a systematic bias in distances recorded with a GPS: the distance between two points recorded with 
a GPS is - on average - bigger than the true distance between these points. This systematic ‘over¬ 
estimation of distance’ becomes relevant if the influence of interpolation error can be neglected, 
which is the case for movement sampled at high frequencies. We provide a mathematical explana¬ 
tion of this phenomenon and we illustrate that it functionally depends on the autocorrelation of 
GPS measurement error (C). We argue that C can be interpreted as a quality measure for move¬ 
ment data recorded with a GPS. If there is strong autocorrelation any two consecutive position 
estimates have very similar error. This error cancels out when average speed, distance or direction 
are calculated along the trajectory. 

Based on our theoretical findings we introduce a novel approach to determine C in real-world 
GPS movement data sampled at high frequencies. We apply our approach to a set of pedestrian and 
a set of car trajectories. We find that the measurement error in the data is strongly spatially and 
temporally autocorrelated and give a quality estimate of the data. Finally, we want to emphasize 
that all our findings are not limited to GPS alone. The systematic bias and all its implications 
are bound to occur in any movement data collected with absolute positioning if interpolation error 
can be neglected. 


1 Introduction 


Global Navigation Satellite Systems (GNSSs), such as the Global Positioning System (GPS), have 
become essential sensors for collecting the movement of objects in geogr aphical space. In movemen t 
ecology, G PS tracking is used to unveil t he migratory paths of birds (Higuchi and Pierr3 . 12005 ), 
elephants ( Douglas-Hamilton et flZ.1 . 12005 1 or ro e deer ( Andrie nko et aLl . l201ll l. In urban studies. 
GPS movement data help detecti ng traffic flows ( Zheng et af.l . l201lh and human activity patterns in 
cities ( Van der Soek e t adl200^ . In transportation research, GPS allows monitor i ng of intelligent 


vehicles ( Zito et aTT 1995 ) and mapping of transportation networks ( Mintsis et al . , 2004[) , to 


name 


but a few application examples. 

Movement recorded with a GPS is commonly stored in form of a trajectory. A trajectory r is 
an ordered sequence of spatio-temporal positions: T=< (Pl,<l),...,(P„,<„) >, with ■■■ tn 
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(|Guting and Schneiden . l2005f) . The tuple {P,t) indicates that the moving object was at a position 
P at time t. In order to represent the continuity of movement, consecu tive positions (Pi,ti) and 


{Pj,tj) along the trajectory are connected by an interpolation function (|Macedo et alJ . l2008f) . 


However, although satellite navigation provides global positioning at an unprecedented accu¬ 
racy, GPS trajectories remain affected by errors. There are two types of error in herent in any kind 
of movement data, measurement error and interpolation error I Schneider . 1999(l . and these errors 
inevitably also affect trajectories recorded with a GPS. 


• Measurement error refers to the impossibility of determining the actual position (P, t) of an 
object due to the limitations of the measurement system. In the case of satellite navigation, 
it reflects the spatial uncertainty associated with each position estimate. 

• Interpolation error refers to the limitations on interpolation representing the actual motion 
between consecutive positions {Pi,ti) and {Pj,tj). This error is influenced by the temporal 
sampling rate at which a GPS records the positions. 


Measurement and interpolation errors cause the movement recorded with a GPS to differ from 
the actual movement of the object. This needs to be taken into account in order to achieve 
meaningful results from GPS data. 

In this article we focus on GPS measurement error in movement data. We show that measure¬ 
ment error causes a systematic overestimation of distance. Distances recorded with a GPS are - 
on average - always bigger than the true distances travelled by a moving object, if the influence 
of interpolation error can be neglected. In practice, this is the case for movement recorded at high 
frequencies. We provide a rigorous mathematical explanation of this phenomenon. Moreover, we 
show that the overestimation of distance is functionally related to the spatio-temporal autocorre¬ 
lation of GPS measurement error. We build on this relationship and provide a novel methodology 
to assess the quality of GPS movement data. Finally we demonstrate our method on two types of 
movement data, namely the trajectories of pedestrians and cars. 

Section introduces relevant work from previously published literature. Section |3] provides a 
mathematical explanation of why GPS measurement error causes a systematic overestimation of 
distance. Section |4] shows how this overestimation can be used to reason about the spatio-temporal 
auto-correlation of measurement error. Section [5] describes the experiment and presents our exper¬ 
imental results, section |6] discusses the results. 


2 Related work 


Since GPS data have become a common component of scientific analyses its quality parameters 
have received considerable attention. These parameters include the accuracy of the estimated 
position, the availability and the update rate of the GPS signal, as well as the continuity, integrity. 


reliability and coverage of the service ( Hofmann-Wellenhof et aJ\ . 2003ll . The accuracy of the 


estimated position (i.e. the expected conformance of a position provided with a GPS to the true 
position, or the anticipated measurement error) is clearly of utmost importance. Measurement 
error and its causes, influencing factors, and scale have been e xtensively dis cussed in published 
literature: measurement error has been shown to vary over time ( Olvnik . I 2 OO 2 II and to be location- 
dependent. Sh adowing effects, for example due to canopy cover, have a significant influence on 


its magnitude ( D’Eon et al. . 200211 . Measurement error is both rando m and caused by ex ternal 


influences, as well as systematic and caused by the system’s limitations (jParent etalj j20 


Measurement error is the result of several influencing factors. According to iLanglev 
these include: 


3). 


(I1997I1 . 
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• Propagation delay: atmospheric variations can affect the speed of the GPS signal and hence 
the time that it takes to reach the receiver; 

• drift in the GPS clock: a drift in the on-board clocks of the different GPS satellites causes 
them to run asynchronously with respect to each other and to a reference clock; 


• Ephemeris error: imprecise satellite data and i ncorrect physica l models affect estimations of 
the true orbital position of each GPS satellite ( Colombo . 1986ll : 


• Hardware error: the GPS receiver, being as fault-prone as any other measurement instrument, 
produces an error when processing the GPS signal; 

• Multipath propagation: infrastructure close to the receiver can reflect the GPS signal and 
thus prolong its travel time from the satellite to the receiver; 


• Satellite geometry: an unfavourable geometric constellation of the satellites reduces the ac¬ 
curacy of positioning results. 

There are several quality measures to describe GPS measurement error, the most common being 
the 95 % radius iR95) , which is de fined as the radius of the smallest circle that encompasses 95 % of 
all position estimates ( Chinl . [r987t) . T he official GPS Performance Analysis Repor t for the Federal 
Aviation Administration issued by the I William J. Hughes Technical Centeii (|20131) states that the 
current set-up of the GPS allows to measure a spatial position with an average R95 of slightly over 
3 meters using the Standard Positioning Service (SPS). The values in the report were, however, 
obtained from reference stations that were equipped with high quality receivers and had unob¬ 
structed views of the sky. It is reasonable to assume that the accuracy would be reduced in other 
recording environments, as measurem ent error depends to a considerable extent on the r e ceiver , 
as well as on the geographic location ( William J. Hughes Technical Center . 2013 : Langlevl 19971. 


This assumption is supported by p ublished literature on G PS accuracy in forests (|Sigrist et al. 


I 999 I I and on urban r oad networks JM odsching et al.. 200dl . as well as on the accuracies of dif¬ 
ferent GPS receivers ( Wing et al\ . 20051 : ^ndbergenl 2009ll . On the other hand, the accuracy of 
GPS can be increased using differential global positioning systems (DGPS) such as the European 
Geostationary Navigation Overlay Service (EGNOS). DGPS esti mate and correct the propagatio n 


delay in the ionosphere, thus yielding higher position accuracies (jHofmann-Wellenhof et al.l . 12003( 1. 


A detailed overview of current GPS accuracy is provided in the quarterly GPS Performance 
Analysis Report for the Federal Aviation Administration. A good introduction to the GPS in 
general, and to its error sourc es and quality parameters in particular, has been provided by 
Hofmann-Wellenhof et al. ( 2003ll . 


The above-mentioned research has mainly focused on describing and understanding GPS mea¬ 
surement error. In addition to this, filtering and smoothing approaches have been proposed for 
recording movement data, in order to reduce the infl uence of erro r s on movement trajecto ries. A 
concise summary of these approaches can be found in IParent et al. ( 2013ll . Jun et al. ( 2006l l tested 
smoothing methods that best preserve travelled distance, speed, and acceleration. The authors 
found that Kalman filtering resulted in the least difference between the true movement and its 
representation. 
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3 GPS measurement error causes a systematic overestima¬ 
tion of distance 

A GPS measurement consists of a spatial component (i.e. latitude (j>, longitude A) and a temporal 
component (i.e. a time stamp t). In this article we will mainly focus on the spatial component. 

The GPS uses the World Geodetic System 1984 (WGS84) as a coordinate reference system. 
For reasons of simplicity it is preferable to transform the GPS measurements to a Cartesian map 
projection, such as the Universal Transversal Mercator (UTM). A transformation from an el- 
l ipsoid (WGS84) to a Cartesian plane (UTM) leads to a distortion of the original trajectories 
( Hofmann-Wellenhof et al.l . l2003ll . For vehicle, pedestrian, or animal movements consecutive posi¬ 
tions along a trajectory are usually sampled in intervals ranging from seconds to minutes. Thus, 
these positions are very close together in space so that the distortion is insignificant for most 
practical applications. Hence, for all the following consideration we assume that the movement is 
recorded in UTM. 

Very generally, a spatial position in UTM is a two-dimensional coordinate 



where x is the metric distance of the position from a reference point in eastern direction and 
y in northern direction. If a moving object is recorded at position P with a GPS, the position 
estimate = {x'^,y^) is affected by measurement error. The relationship between the true 
position and its estimate is very trivially 


P^=P + ep, (2) 

where £p is the horizontal measurement error expressed as a vector in the horizontal plane. £p 
is drawn from £p, the distribution of measurement error at P. We adopt the convention used 


"ft K 

by iGodlinsf et all ([20080 to denote random variables with upper case letters and their numerical 
values with lower case letters. 

We now provide a detailed mathematical explanation of why measurement error causes a sys¬ 
tematic overestimation of distance in trajectories, if interpolation error can be neglected. Figure [T] 
illustrates the problem statement in a simplified form. Gonsider a moving object equipped with 
a GPS. The moving object travels between two arbitrary positions P and Q. Let do = d{P, Q) 
denote the Euclidean distance between these, henceforth referred to as reference distance. The 
object always moves along a straight line, consequently interpolation error can be neglected. The 
movement of the object can be described by the following five steps, which correspond to the 
subplots in Figured] 


1. The moving object starts at P. The GPS obtains the position estimate P'^ with measure¬ 
ment error ep, which is drawn from £p. 

2. The moving object travels to Q. The GPS obtains the position estimate Q'^ with measure¬ 
ment error sq, which is drawn from £q. The distance between the two position estimates is 
calculated: d™ = d{P^,Q^). 

3. The moving object returns to P. The GPS obtains a position estimate and a new d™ is 
calculated. 

4. Steps 2 and 3 are repeated n times, where n is an infinitely large number. 

5. After n repetitions, the position estimates scatter around P and Q with measurement error 
£p and £q- 
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Figure 1: A moving object equipped with a GPS travels between two arbitrary positions 


We claim that measurement error propagates to the expected measured distance and to 

the expected squared measured distance E((i™) between the position estimates. More specifically, 
measurement error yields E(ci'") > do as well as E((i™) > dp. 

We are now going to rigorously prove this claim. To do so, we simplify notation, write £p — 
{Xi,Yi) as well as £q = {X 2 ,Y 2 ), and assume that there is no systematic bias, i.e. we have 
E(Xi) = E(A 2 ) = E(li) = E(y 2 ) = 0. Since neither translations nor rotations affect distances 
between points we may, without loss of generality, consider P = (0,0) and Q — ((io,0). Since 
linear transformations (like rotations) preserve expectation, rotating errors with expectation zero 
results in errors having expectation zero too. Having this we can now formulate the following first 
re sult for t he exp ected squared distance E(d^(P™', Q™')). For mathematical background we refer 
to [m en ki ll201.4 . Notice that no assumptions (like absolute continuity or normality) about the 
underlying error distributions are needed, i.e. the result holds in full generality. 


Theorem 1 Suppose that do > 0, that P = (0,0), and that Q = (dojO). Let Xi,X 2 both have 
distribution function F and variance u\, and Yi , Y 2 both have distribution function G and variance 
a\. Furthermore assume that E(Xi) = E(A 2 ) = E(Yi) = E(y 2 ) = 0. Then the following two 
conditions are equivalent: 

1. E(d^) = E(d2(P"»,Q’")) > dl 


2 . min{C'ou(Ai, A2), C'ou(Fi, 1^2)} < 1 

In other words: the expected squared distance E(d™) is strictly greater than d§ unless the errors 
fulfil Xi = X 2 and Yi = Y 2 with probability one (which describes the situation of always having 
identical errors in P and Q). 


Proof: Calculating E(d^(P™', Q™')) and using the fact that C'ou(Ai, A 2 ) < and C'o?;(Yi, ^ 2 ) < 
CTy directly yields 


E(d2(P"^,Q"*)) = 


> 


E(do + A2 - X{f + E(Y 2 - Yif 

dl + E(A2 - X^f + E(Y2 - Yif = dg + Var{X 2 

dl + + 2ct^ - 2C'ou(A i,A 2) - 2C'ou(Yi,y2) 


Ai) + Var{Y2 


Yi) 

( 3 ) 
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Having this it follows immediately that E(d^(P"^, Q™')) = ci§ if and only if Cov{Xi, X 2 ) = cf\ 
and Cov{Yi,Y 2 ) = ay, which in turn is equivalent to the fact that Xi = X 2 and Yi = Y 2 holds 
with probability one. ■ 

In general one is, however, interested in the expected distance E{dX) := E(d(P™', and 

not in the expected squared distance. Since, in general, E(Z^) > dp need not imply E(|Z|) > dp 
for arbitrary random variables Z, a different method is used to prove the following main result: 

Theorem 2 Suppose that the assumptions of Theorem]^ hold. Then the following two conditions 
are equivalent: 

1. E(d’”) = E(d(P"*,Q”")) > do 

max{P(ri ^ y 2 ),P (^2 -^1 < -do)} > 0 

In other words: the expected distance E(d"*) is strictly greater than the true distance dp unless the 
errors fulfil Yi = y 2 with probability one and f’{X 2 — Xi < —do) = 0 holds. 

Proof: Obviously we have 

x/(do + Y 2 - Xi)2 + (Y 2 - Yi)2 > |do + Y 2 - Xil (4) 

Setting Z := X 2 — Xi implies E(Z) = 0. Assume now that P(Y < -dp) > 0 holds. Then the 
desired inequality follows immediately from 

E|Z + do| = [\z + do\dF^=[ (z + do)dP^+ / -(z + do)dP^ 

«/]R J[—do,oo] J{ — oo, — do) 

= [ iz + dp) dP^ + (-2) [ {z + do) dP^ (5) 

d]R d( —00, —do) 

^---' '-V-' 

=do >0 

> do. 


In case we have F(Z < —do) = 0 but P(Yi ^ Y 2 ) > 0 holds, then Inequality |4] is strict with positive 
probability so we get 


E(d”^) = E(V(do + A 2 - Xi)2 + (Y 2 - Yi)2) > E(|do + A 2 - All) = E(|Z + dpi) = do. 

Altogether this shows that the second condition of Theorem [5] implies the first one. 

To prove the reverse implication, assume that max{P(Yi ^ y 2 ),P(A 2 — Ai < —do)} = 0. Then, 
firstly, the left and the right hand-side of Inequality |4] coincide with probability one, so E(d'") = 
E(\/(do -I- A 2 — Ai)2 -I- (Y 2 — Pi)^) = E(Mo + -^2 — All) holds. And secondly, directly applying 
Equality [ 5 ] yields E(|Z -|- do|) = dp, which finally shows E(d"‘) = dp. ■ 

Remark 1 It is worth mentioning that Theorem [5] has several interesting (and partially surpris¬ 
ing) consequences: Whenever the errors in a;-direction are unbounded (like in the case of normal 
distributions) the expected distance is always strictly greater than the true distance do- The 
same holds whenever the errors Yi and Y 2 in j/-direction do not always coincide - a very realistic 
assumption for GPS trajectories. 


We want to underline that Theorem [T] and H] hold in full generality for arbitrary distributions of 
GPS measurement error. Although GPS measurement error is often assume d to have a bivariate 
normal distribution a n d to be ind ependent in both the x- and y-direction (|Jerde and Visscher . 
2005 ; Bos et al . , 20081) , IChinI ( 19871) puts forward convincing arguments why this is very likely not 


the case. Hence, the general validity of our findings is relevant. 
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For reasons of simplicity, we have assumed that £p and Sq follow the same distribution function 
and that there is no systematic bias, i.e. £p is centred around P and Sq around Q. This 
assumption is generally acknowledged for in the literature. It b uilds, for example, the basis for 
algorithms to extract road maps from GPS tracking data (e.g. Wang et al\ . l2014ll . Roads are 
assumed to be located where the density of the GPS measurements is the highest. Also Figure |4] 
shows that this assumption is indeed realistic for real-world GPS data. However, even a systematic 
bias does not necessary restrict the validity of our argument. Let us assume that ]E(Ai) = E(A 2 ) 7 ^ 
0 and E(Yi) = E(y 2 ) ^ 0, i.e. the mean of the error distribution has shifted away from P and Q 
respectively. As the shift is the same for £p and 5 q, the influence on distance calculations cancels 
out. Theorem [T] and [5] still hold. The validity of our proof is restricted only if E(Ai) 7 ^ E(A 2 ) or 
E(Fi) 7 ^ E(y 2 )- This implies that the mean of the error distribution changes abruptly between P 
and Q. As - in practice - P and Q are very close in space, this scenario is not realistic for GPS 
measurement error. 


4 How big is the overestimation of distance and why is this 
relevant? 

In the previous Section we proved that distances recorded with a GPS are on average bigger than 
the distances travelled by a moving object, if interpolation error can be neglected. In this Section 
we provide an equation for OD, the expected overestimation of distance. Moreover, we identify 
three parameters that influence the magnitude of OD. First, let us define OD with the help of 
Equation [31 

OD=E{d^)^ -do = {dl + -h 2 al- - 2 Goz;(Ai, A2) - 2 Cov{Yi,Y 2 ))i - do. 

From this follows that OD is a function of three parameters: 

1. do, the reference distance between P and Q 

2. Vavgps = 2aj^ + 2CTy, a term for the variance of GPS measurement error 

3. C = 2 Cov(Xi , X 2 ) — 2Cov{Yi , 12)1 a term for the auto-correlation of GPS measurement error. 
C expresses the similarity of any two consecutive position estimates. If C is big, consecutive 
position estimates are affected by similar GPS measurement error (see also Figure 121). 

We can now simplify notation and write 

OD = {dl + Vavgps -C)^ - do. (6) 

The influence of the three parameters on OD is further illustrated in Figure |21 OD is small if 
the reference distance is big, the variance of GPS measurement error is small and the error has 
high positive autocorrelation. OD is big if the reference distance is small, the variance of GPS 
measurement error is big and the error has high negative autocorrelation. 

To understand the magnitude of OD in real-world GPS data, let us assume for a moment that 
there is no autocorrelation of GPS measurement error, i.e. (7 = 0. Moreover, let us assume that 
the variance of error is the same in x- and y- directions, i.e. cr^ = = Oy and Vaxgps = dcr^. We 

can now visualize the relationship between OD, do and a. Figure [31 a shows that OD increases as 
the spread of GPS measurement error (cr) increases; dp is assumed to be constant. For a constant 
do of 5m, for example, and a = 2m the overestimation of distance roughly equals 2m (yellow line). 
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Figure 2: Overestimation of distance {OD) and its influencing parameters. 


When a increases to 4m the overestimation of distance increases to 4m. Figure [3] b shows that OD 
decreases as do increases; a is assumed to be constant. For a constant a of 3m, for example, and 
do = 5m the overestimation of distance equals around 3m (black line). When do increases to 10m 
the overestimation of distance decreases to 2m. 

Remember that Figure |3] shows the influence of Vavgps if there is no a utocorrelation of GPS mea¬ 


surem e nt error. Thi s is no t ver y realistic for real wor ld GPS data. In fact. lEl-Rabbanv and Kleusbere 
(120031) : I Wang et al\ (120021) and iHowind et all ()l999l) show that GPS is affected by both spatial and 
temporal autocorrelation. This means that position estimates taken close in space and in time 
tend to have similar error. 

How big is the autocorrelation of GPS measurement error? Let us reformulate Equation |6] and 
solve for C: 


K{dT) 

C = do — {OD + do)'^+Vargps. (7) 

This implies that we can calculate the autocorrelation of GPS measurement error if OD, Vavgpg 
and do are known. Things become interesting if we consider what autocorrelation really means in 
the context of GPS positioning. In Figure [5J in the bottom left cell, the position estimates 
and Q™' are highly auto-correlated and, hence, very similar. This leads to the effect that d'" is 
very similar to dp. In fact, this applies not only to distance, but to other movement parameters as 
well. Direction, speed, acceleration or turning angle must all be similar to the ‘true’ movement of 
the object if they are derived from highly auto-correlated GPS position estimates. Gonsequently, 
C describes how well a GPS captures the movement of an object, if interpolation error can be 
neglected. Or in other words, C is a quality measure for GPS movement data. 


































a 


b 




Figure 3: The overestimation of distance (OD) increases as the spread of GPS measurement error 

(a) increases, the reference distance (do) is constant (a); OD decreases as do increases a is constant 

(b) . 


5 Assessing the quality of GPS movement data 

Real world GPS data are affected by spatial as well as temporal autocorrelation (lEl-R,abbanv and Kleusb^ . 


2003 : Wang et al . 20021 Howind et al . lOOOl l. Spatial autocorrelation implies that GPS measure¬ 


ment error is not independent of space. Measurements obtained at similar locations will have 
similar error. Temporal autocorrelation implies that GPS measurement error is not independent 
of time. Measurements obtained at similar times will have a similar e rror due to similar atmo¬ 
spheric conditions and a similar satellite constellation ( Bos et al. . 2008h . We carried out a simple 
experiment to visualize temporal autocorrelation in real-world GPS data. We placed a GPS unit 
at a known position P and recorded about 720 position estimates over a period of about six hours 
at a sampling rate of 1/30 Hz. The resulting distribution is centred around P with an R95 of 
about 3 m (Figured a). If only those position estimates are displayed that were recorded within a 
certain time interval, GPS measurement error reveals itself to be highly auto-correlated. Figure |4] 
b, for example, shows only those position estimates that were obtained within periods covering 
5 minutes before and after 

In this Section we build on the relationship described in Equation [3 and show the spatial and 
temporal autocorrelation in two sets of real-world GPS movement data. In the first experiment we 
identified to what degree a set of pedestrian movement data was affected by spatial and temporal 
autocorrelation. In the second experiment we derived the spatial autocorrelation in a set of car 
movement data. Based on this we tried to assess how well the GPS captured the movement of the 
car. 


5.1 Experiment 1: Pedestrian trajectories 

Experimental setup 

For the first experiment, we equipped a pedestrian with a GPS. The pedestrian walked along a 
reference course with a well-established reference distances do- The movement of the pedestrian 
was recorded with a QSTARZ:BT-Q1000X GPS unil[3 with ‘Assisted GPS’ activated. 

^for specifications, please refer to: http://www.qstarz.eom/Products/GPS/20Products/BT-Q1000.html 


9 














































Figure 4: The distribution of GPS measurement error at position P (a). Revealing the temporal 
autocorrelation of GPS measurement error (b). The movement of a pedestrian around a reference 
course (c). 


Rather than using a high-quality GPS we collected all data with a low-budget GPS, a type of 
GPS common for recording movement data. We deliberately treated the GPS as a ‘black box’. 
This implies that the algorithm to calculate the position estimates from the raw GPS signal was 
not known. Moreover, we considered that it was sufficient to use only a single GPS unit, as the 
aim of the experiment was not to investigate the quality of the particular GPS, but to show the 
usefulness of our approach. 

The reference course was located in an empty parking lot to avoid shadowing and multi-path 
effects. We staked out a square with sides that were 10 m long and had markers at one meter 
intervals. A square was used in order to allow distance measurements to be collected in all four 
cardinal directions (approximately). The distance between the markers was used as a reference 
distance do. 

The GPS position estimates were obtained by walking to the reference markers in turn and 
recording the position, moving around the square until all positions of the markers had been 
recorded. Position estimates were only taken at the reference markers, and only when the recording 
button was pushed manually. Two consecutive position estimates were taken within three to five 
seconds. A full circuit around the square took approximately between two and three minutes and 
resulted in 40 positions being recorded. A total of 25 circuits around the square were completed, 
without any breaks. This resulted in 1000 GPS positions being collected in approximately one 
hour. 

In pre-processing distance measurements d'" were calculated between the position estimates 
and later compared to do the reference distance between the markers. Then the average measured 
distance d™ was calculated and from this OD = d'" — do and (7 = dg — d™ -I- Vargps were derived. 
OD and C are estimators for OD and C. 

We decided not to derive ax and ay from observational data, but to set ax = ay = 3m. 
Hence, Vargps is not the observed variance of GPS measurement error, but a reference value to 
which OD is later compared to. Consequently, our results do not show the exact value of C, but 
provide an estimate of C with respect to Vargps. 

We increased the spatial separation between two position estimates of the pedestrian to illus¬ 
trate the influence of spatial autocorrelation. Then we increased the temporal separation between 
two position estimates to illustrate the influence of temporal autocorrelation. 
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Figure 5: Overestimation of distance {OD) and spatial autocorrelation of GPS measurement error 
(C) in the pedestrian movement data. 


Results 

In contrast to the theoretical findings in Figure |3l overestimation of distance tended to increase as 
the reference distance dg increased. This was due to a decrease in the spatial autocorrelation of 
GPS measurement error. With increasing spatial separation of the position estimates, measurement 
error became less auto-correlated. Figure [S] shows the relationship between the reference distance 
do and OD (black dots) as well as C (black crosses). 

We wanted to illustrate that the overestimation of distance was not caused by a small number 
of extreme outliers. Figure [6] shows the histogram of d™ — dg for dg = 1 m (a), and for dg = 5 m (b) 
and their fit to a Gaussian distribution. Both histograms follow a Gaussian distribution cr^) 

rather well and outliers are almost non-existent. Note that pd and in Figure [6] refer to the 
values of the fitted Gaussian distribution and not to the empirically derived frequency. 

In order to illustrate the temporal autocorrelation in GPS measurement error, we calculated 
the distance between non-consecutive position estimates around the square. One example is the 
distance between two position estimates, where the second one was obtained one circuit after the 
first later. The reference distance between the markers remained the same, e.g. dg = 1 m, but the 
position estimates were recorded within a longer time interval At. Figure [7| shows the relationship 
between At and OD (black dots) as well as C (black crosses) for a reference distance dg = 1 m. 
OD increase with longer time intervals. The sharpest increase occurs between measurements that 
were taken promptly and those taken after about 2^ minutes. After 40 minutes the curve levels 
out. This increase of OD was caused by the temporal auto-correlation of measurement error. 
For measurements taken within several seconds, measurement error appears to be strongly auto- 
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Figure 6: Histogram of the difference between measured and reference distance (d™ — dp) for 
dp = 1 m (a) and dp = 5 m (b) 


correlated. However, anto-correlation falls sharply for measurements taken within 2^ minutes. 
From then on C gradually decreases as At increases; again the curve levels out at about 40 
minutes. 

The data for the above experiment were calculated with a GPS for which the algorithm to 
calculate the position estimates from the raw GPS signal was not known. This raises the legitimate 
question, whether the results were produced by a smoothing algorithm rather than the behavionr 
of the GPS. Let us assume that the GPS used a smoothing algorithm. In simplified form, the 
current position estimate is then calculated from the last position estimate, the current GPS 
measurement and a movement model. For movement with constant speed and direction, smoothing 
yields trajectories that represent the true movement very accurately. However, sudden changes in 
movement, i.e. a sharp turn, are not followed by the trajectory. The current measurement implies a 
sharp turn, however, the movement model does not. Thus, the sharp turn becomes more elongated, 
the overestimation of distance increases. However, we did not find any support for an increase in 
the overestimation of distance after a sharp turn. This can also be seen in Figure |4]b. 

5.2 Experiment 2: Car trajectories 

In the first experiment the reference distance cLq was staked out along a reference course. For 
obvious reasons this is not possible for recording the movement of a car. Hence we derived dp from 
speed measurements recorded with a car’s controller area network bus (CAN bus). 

Experimental Setup 

We equipped a car with a GPS unit and tracked its movement for about 6 days. The car moved 
mostly in an nrban road network at rather low speeds (average: 25 km/h). The temporal sampling 
rate of recording was 1 Hz. For the CAN bns measurements, a sensor recorded the rotation of 
the car’s drive axle, from which dp was inferred. Thus dp is the distance travelled by the car 
according to the CAN bus. For the same phases of movement we compared dp to d™, the distance 
travelled by the car according to the GPS position estimates. As in the hrst experiment, we set 
o'x = <Jy = 3m and calculated Vargas. 
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Figure 7: Overestimation of distance OD and temporal autocorrelation of GPS measurement error 
(C) in the pedestrian movement data. 


The data were first pre-processed and cleaned. Parts with very high speed (above 140 km/h) 
and very rapid acceleration (above 5 m/s^) were removed. Although the data consisted mostly of 
the car’s forward movements, there were also periods when it was either stationary or reversing in 
a parking lot. The data may also have included some periods during which shadowing caused a 
loss of the GPS signal (for example when driving in a tunnel). We therefore applied a simple mode 
detection algorithm to remove any such periods. The algorithm evaluates speed and acceleration 
along the trajectory and distinguishes segments that m ost probably r eflect driving behaviour from 
those that are likely to reflect non-driving behaviour ( Zheng et '^■ l2010ll . Using the algorithm 
we were able to include only long phases of continuous driving, sampled at a continuous sampling 
frequency of 1 Hz. Following this pre-processing a total of about 195 km of car trajectories remained 
for analysis. 


Results 

Figure [8] shows that the autocorrelation of GPS measurement error decreased as the spatial sepa¬ 
ration between two consecutive position estimates increased. Nevertheless, C in Figure [8] is always 
positive. This can be interpreted as a quality measure for the movement data. Consecutive position 
estimates were affected by less variance than initially suggested by Vargpg. 

Although the results in Figure [5] are similar to those obtained from the pedestrian movement 
data, they contain outliers. We believe that these outliers occur due to two reasons. First, the 
data comprise relatively few distance measurements for big do because of the generally low speed 
of the car. Second, we could not guarantee a full temporal synchronization of both measurement 
systems (GPS and CAN bus). In other words, do and might relate to slightly different time 
intervals. We found this lag to be around one second. We believe that this insight is important 
for the practical application of Equation 0 In order to provide valid results it requires both a 
signihcant number of distance measurements as well as a proper synchronisation of reference and 
measured distance. 


13 














Figure 8: Overestimation of distance {OD) and spatial autocorrelation of GPS measurement error 
(C) in the car movement data. 

6 Discussion and Outlook 


In this paper we identified a systematic bias in GPS movement data. If interpolation error can 
be neglected GPS trajectories systematically overestimate distances travelled by a moving object. 
This overestim ation of distance has previously been noted in the trajectories of fishing vessels 
( Palmer! . I 2 OO 8 I I . For high sampling rates the distance travelled by the vessel was overestimated due 
to measurement error, while for lower sampling rates it was underestimated due to the influence of 
interpolation error. We provided a mathematical explanation for this phenomenon and showed that 
it functionally depends on three parameters, of which one is C, the spatio-temporal autocorrelation 
of GPS measurement error. We built on this relationship and introduced a novel approach to 
estimate C in real-world GPS movement data. In this Section we want to discuss our findings and 
show their implications for movement analysis and beyond. 

In the era of big data, more and more movement data are recorded at hner and hner intervals. 
For movement recorded at very high frequencies (e.g. IHz) interpolation error can usually be 
neglected. Hence OD is bound to occur in these data. However, this does not mean that high 
frequency movement data are of low quality, quite the opposite is true. Using the relationship 
between C and OD we showed experimentally that GPS measurement error in real world trajecto¬ 
ries was affected by strong spatial and temporal autocorrelation. In other words, if the data were 
recorded close in space and time they captured the movement of the object better than if they 
were further apart. 

Autocorrelation is important for movement analysis in many aspects. An appropriate sampling 
strategy for recording movement data, for example, should consider the influence of measurement 
error and address spatial and temporal autocorrelation. Since autocorrelation can be interpreted as 
a quality measure, it allows to reveal the performance of different GPS receivers in different record¬ 
i ng en vironments. Moreover, autocorrelation has implications for simulation. Laube and Pur\^ 
(|201l[l performed a simulation to reveal the complex interaction between measurement error and 
interpolation error and their effects on recording speed, turning angle and sinuosity. Their Monte 
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Carlo simulation assumed GPS errors to scatter entirely randomly between each two consecutive 
positions. Our approach allows to verify whether this assumption is realistic. 


Where to find a reference distance? 


For practical applications the biggest limitation of our experiments is their dependency on a valid 
reference distance. The moving object must traverse the reference distance along a straight line 
and without interpolation error, and at a precisely known time. Moreover, a large number of 
measurements has to be collected, since C is derived from the expectation value of a random 
variable. 

This limitation leads to a possibly interesting application of our findings, where the reference 
distance is derived from the GPS point speed measurements. Point speed measurements are cal¬ 
culated from the instantaneous deriva tive of the GPS signal using the Doppler effect. Point speed 
is very accurate I Bruton et all . Il999ll and usually part of a GPS position estimate. Hence, for 
high sampling rates (e.g.l Hz) point speed measurements can be used to infer the distance that 
a moving object has travelled between two position estimates. This distance is not affected by 
the overestimation of distance effect and could serve as a reference distance. Thus, GPS could be 
compared to itself to reveal the spatio-temporal autocorrelation of the position estimates. This 
approach would not require any other ground truth data. However, its feasibility and usefulness 
are yet to be tested. 

Finally, we want to underline that our findings are not only relevant for GPS. The overestimation 
of distance is bound to occur in any type of movement data where distances are deduced from 
imprecise position estimates, of course only if interpolation error can be neglected. 
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